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Abstract 

A simple uniform approximation of the logarithmic derivative of the ground state 
eigenfunction for both the quantum-mechanical anharmonic oscillator and the double- 
well potential given by V = m 2 x 2 + gx A at arbitrary g > for m 2 > and m 2 < 0, 
respectively, is presented. It is shown that if this approximation is taken as unper- 
turbed problem it leads to an extremely fast convergent perturbation theory. 
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For the last fifty years the quantum one-dimensional anharmonic oscillator (m? > 0) as well as 
the double-well potential problem (m 2 < 0) described by the Hamiltonian 

d 2 

H = + m 2 x 2 + gx i , (1) 

permanently attracted a lot of attention. The interest to these problems ranges from various 
branches of physics to chemistry and biology. It can not be an exaggeration to say that after 

n 

the seminal papers by C. Bender-T.T. Wu at 1969-1973 d] in near thousand of physics articles 

the problem (^Q) was touched in one way or another. This seemingly simple problem revealed 

extremely rich internal structure which looks intrinsic for any non-trivial problem of quantum 

mechanics and even for quantum field theory. In particular, the obvious failure of the perturbation 

theory in powers of g due to its asymptotic (divergent) nature in Q pushed the development of non- 

perturbative methods. In practice, the anharmonic oscillator (^Q) served always as a test-ground for 

non-perturbative methods. Another important feature of the anharmonic oscillator Q is related 

to the fact that it can be interpreted as one-dimensional quantum field theory gcp 4 with zero spacial 

dimension, (1,0). Hence, a study of Q can be insightful to a realistic four-dimensional quantum 

field theory gcj) 4 . The present article is devoted to the construction of a uniform approximation of 

the ground state eigenfunction of (Q) in x-space which would continue to be valid for different values 

of the coupling constant g and the parameter m 2 . If such an approximation is constructed the 

energy of the ground state as well as different average values can be calculated with a guaranteed 

accuracy. We are not aware of previous attempts to construct a uniform approximation of the 

eigenf unctions. 

Take the Schroedinger equation for Q 

d 2 ^ f +oc 
--j~2 + m2x2 ^ + 9 X ^ = E(m 2 ,g)^ , J \^>\ 2 dx < oo . (2) 

For m 2 > the potential in © is a single-well potential, which describes the celebrated quartic 
anharmonic oscillator. For m 2 < the potential in © is a two-well potential. It describes 
another celebrated potential - the so-called double- well potential (also known as the Higgs potential, 
Lifschitz potential). It is easy to check that for the energy in © the Symanzik scaling relation 
holds 

E(m 2 ,g)=g^E^l) , ^;m 2 , 5 ) = *(V /6 ; l) ■ 
This manifests that the original problem © is essentially a single-parametric problem. 
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Eigenfunctions of (J2J) are sharply changing functions in x being characterized by a power-like 
behavior at \x\ — > and an exponentially-decaying one at |x| —* oo. In order to make them smooth 
we introduce the exponential representation for eigenfunctions 

^j(x) = e -v(x) _ ( 3 ) 

Following the oscillation (Sturm) theorem the function (p should have logarithmic singularities 
at real x which correspond to nodes of the wavefunction. Its regular part which appears after 
subtraction of those singularities should be slow-changing function of the power-like behavior at 
both small and large distances. After substitution of (jHJ) into we get a Riccati equation 

y'-y 2 = E- m 2 x 2 - gx 4 , y = cp' = (log*(a?))' , y(0) = . (4) 

This famous equation serves as a basis for developing the WKB approximation scheme. Due to 
symmetric nature of the r.h.s. of Q) x — > —x the function y is odd. Its structure for the nth 
excited state should be 

n 1 

y = + y^ x ) . ( 5 ) 

i=l 

where xi is the position of ith node. The regular part y reg {x) should be a non-singular at real x 

odd function which vanish at x = Xi, i = 1, 2, . . . n. It is evident that all non-zero Xi come in pairs: 

once we have a node at x = x% always there exists a node at x = — x,. 

Now we consider the case of the ground state. According to above, in the phase © n = and 

singular part of y is absent, y = y reg . Hence, the function y has no singularities at real x. It is 

easy to find its asymptotic behavior 

1/2 i \x\ 1 4gE + m 4 1 m 2 1 

V = 9 Xlxl + ^- + x-^gJJ^x]x]-^ + --- at N^°°' ( 6 ) 

while 

E 2 -m 2 o 2E(E 2 -m 2 )-3g «. 

y = Ex^ x A H -x 5 + ... atxUO. (7) 

y 3 15 11 w 

One can demonstrate that for any excited state the regular part y reg has similar formal expansions 

at \x\ — > 0, oo. 

Let us develop a certain iterative procedure (perturbation theory) for solving the Riccati equa- 
tion @. From the point of finding the wave function it will be a multiplicative perturbation theory 
unlike a standard additive perturbation theory. Such a perturbation theory was developed for the 

3 



first time by Price 



and then it was rediscovered many times. Eventually, it was called the 
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'Logarithmic Perturbation Theory' (for the history remarks and discussion see [3( and references 
therein). 

As a first step we choose some square-integrable function and calculate its logarithmic 
derivative 

y = (log*o)' = 5 • (8) 
It is clear that is the exact eigenfunction of the Schroedinger operator with a potential 

Vo = ^ = y 2 -y' , (9) 

where without a loss of generality we put their eigenvalue equals to zero, Eq = 0. It is nothing but 
a choice of the reference point for eigenvalues. Now we can construct a perturbation theory for 
Riccati equation taking an d yo, Vq as zero approximation, which characterizes the unperturbed 
problem. One can write the original potential V = m 2 x 2 + gx A as a sum, 

V = V + (V - Vq) = V + V l: (10) 

thus, taking a deviation of the original potential from the potential of the zero approximation as a 
perturbation. We always can insert a formal parameter A in front of V\ and develop a perturbation 
theory in powers of A, 

oo oo 

E = J2^Ek , y = ^A fc y fc , (11) 

k=0 k=0 

putting A = 1 afterwards. Perhaps, it is worth emphasizing that in spite of the fact that we study 
iteratively the equation <|l)l. in general, this perturbation series has nothing to do with a standard 
WKB expansion. By substituting pip into @ we arrive at the equations which defines iteratively 
the corrections 

y'k - 2 yoyk = E k -Q k , (12) 

where 

Qi = V 1 , 

k-l 



Qk = -/jVi- Vk-i , k = 2,3,. 



i=l 



It is interesting that the operator in the l.h.s. of p2p does not depend on k, while Qk in the r.h.s. 
can be interpreted as a perturbation on the level k. The solution of p2p can be found explicitly 
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and is given by 



(13) 



Vk = 




— oo 



(E k -Q k )*l dx 



(14) 



It is easy to demonstrate that 



2/1 1 < Const 



(15) 



provides a sufficient condition for this perturbation theory ((TTj) to be convergent |3j. Note that 
this condition is very rough and can be strengthened. 

The first two terms in the expansion of energy 1)11(1 in the above-described perturbation theory 
admit an interpretation in the framework of the variational calculus. Let us assume that our 
variational trial function ^o(x) is normalized to 1. We can calculate the potential Vq where *&o(x) 
is the ground state eigenfunction and even put Eq = (see a discussion above). Pure formally, we 
construct the Hamiltonian Hq = p 2 + Vq for which Hq^q(x) = 0. The variational energy is equal 
to 



Of course, &o(x) could depend on free parameters. In this case both Vo and V\ depend on param- 
eters as well. Minimization of E var with respect to the parameters can be performed and the vari- 
ational principle guarantees that E var gives upper bound to the ground state energy. This simple 
interpretation ()16j) reveals a fundamental difference between perturbation theory and variational 
calculus. Variational estimates can be obtained independently on the fact that the perturbation 
theory associated with *$>o(x) is convergent or divergent. However, it seems natural to remove this 
difference by requiring a convergence of the perturbation series. In this case by calculating the 
next terms E2, -E3 , . . . in one can estimate the accuracy of variational calculation from one side 
and improve it iteratively from another side. An immediate criteria how to choose ^>o(x) in order 
to get a convergent perturbation theory is to have the perturbation potential V\ to be subordinate 
with respect to the non- vanishing potential of zero approximation Vq, 



E, 



'var 




E + E 1 (V 1 = V-V )>E exact . 



(16) 




for \x\ > R 



(17) 
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This requirement has extremely non-trivial physical implication: in order to guarantee a conver- 
gence of perturbation theory a domain where the wavefunction is exponentially small should be 
reproduced as precise as possible. The same time a description of a domain where the wavefunction 
is of the order 1 is not important. It contradicts to a straightforward physics intuition and under- 
lying idea of variational calculus which, in particular, requires a precise description of the domain 
where the wavefunction is of the order 1. Needless to say that namely the latter domain gives a 
dominant contribution to the integrals which define the energy in the variational calculations. 

In our approach the main object for study is the logarithmic derivative of the wavefunction y 
(see ©). Since y is antisymmetric, we will construct different interpolations of y between x = 
and x = oo (see ©,©)• In order to fulfil the requirement of the convergence Q17j) it is enough to 
take into account in these interpolations the leading term of the asymptotics at \x\ — > oo. Then 
the interpolation we have built is taken as zero approximation in our perturbation theory 1)11(1 . 

1. The simplest interpolation can be written as follows 

yo = ax + b^/gx\x\ , (18) 

where a, b are parameters. These parameters can be fixed either by taking them as variational and 
making a minimization of (fl6|) . a m i n ,b m i n , or following the idea to reproduce exactly the leading 
asymptotic behavior of y at x = oo (see (©), which requires to put b = 1. 
The ground state eigenfunction which corresponds to (|18(1 

4 1} =exp{-^-&f W 3 } , (19) 
is the exact one for the potential 

V = a 2 x 2 + b 2 gx* - 2b^/g\x\ (1 - ax 2 ) , E = a . 
The perturbation potential V\ = V — Vq is of the form 

Vx = (m 2 - a 2 )x 2 + (1 - b 2 )gx^ + 2b^j\x\{\ - ax 2 ) , 

where the first two terms would disappear if we place a = ±m, and 6 = 1. In this case the first 

two terms in the perturbation theory are 

m \™ x(\-mx 2 )e- mx2 -^dx , s 

J5« =E +E 1 =m + 2^g ^ — . (20) 

This expression leads to sufficiently high relative accuracy < 10~ 2 (comparing to the accurate 
numerical results) for any g > and value of m 2 . It can be confirmed by calculation of the second 
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correction E<i- However, a slight modification of of the interpolation (|18[) by including a term 1/x, 
which is presented in Q, in a form 



where c is a variational parameter, immediately leads to a drastic increase in accuracy (see Table 1). 
According to a quite simple, straightforward analysis of the second correction E2 the variational 
energy deviates from exact one in < 10~ 3 in relative units for both anharmonic oscillator and 
double- well potential for studied values of the parameters in (1) (see Table 1). When the second 
correction is taken into account the relative deviation reduces in two orders of magnitude becoming 
< 10 -5 . It is worth mentioning that the quality of the approximation is reflected in the fact that 
bmin deviates from the exact value b = 1 in several percent. 

The ground state function is symmetric w.r.t. x — > —x. Hence, it has an extremum at x = 0. 
For any fixed m 2 < there exists a value g cr it such that for g > g cr u this extremum is a maximum, 
otherwise a minimum. It is easy to find out that the critical point g = g cr u corresponds to the 
vanishing ground state energy, E = 0. Using the function 1)21(1 it was calculated the critical value 
gcrit = 0.302405 for m? = —1. It is quite interesting from physical point of view that for a family of 
double-well potentials with fixed g there exists a domain > m 2 > (m 2 ) cr it where the ground-state 
eigenfunction has the maximum at the origin, which corresponds to the position of the unstable 
equilibrium similar to what takes place for the single- well case. For example, if g = 1, the value of 
(m 2 ) cr it = —2.219597. It implies that the particle in such a potential with the ground state energy 
above the barrier, E > 0, somehow does not feel the existence of two minima. 

2. The expression E^ (|20[) reproduces the harmonic oscillator energy in the limit g — ► 0. At 
g = the energy E^ l > has a singularity as it must be. However, the expansion of (|20j) in powers of 
g contains besides the integer powers in g also half-integer powers which must not be present in the 
formal expansion of the energy in g (see (|11|) . Also the interpolation (|18|) reproduces the leading 
terms only in the expansions at |x| — > 0, 00. It is a definite drawback of the interpolation (|18|) as 
well as (|21|) and it should be fixed. One of the simplest ways to fix it is to take the function 




(21) 




(0,0) 




(22) 



where A, a, b, c, d are variational parameters. If the parameters are chosen to be 



b 



4 



3 



a 



d 2 2 

= - + m 



(23) 
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= 2 


9 = 


= 1 


9 = 


= 2 




2 

m 


= 1 


2 

m 


= 


2 

m = 


= -1 




6 = 1 


bmin 


6 = 1 


bmin 


b = 1 


bmin 






=1.09320 




=1.13049 




= 1.19610 




1.6076526 


1.6076150 


1.0605130 


1.0604314 


1.0299142 


1.0296682 


E2 


-0.0001113 


-0.0000735 


-0.0001513 


-0.0000692 


-0.0003560 


-0.0001070 


£(2) 


1.6075413 


1.6075415 


1.0603617 


1.0603622 


1.0295581 


1.0295612 


Ref.[4] 


1.60754130 


1.06036209 


1.02956085 



TABLE I: The Hamiltonian (1): Variational energy of the ground state and its corrections in 
perturbation theory (|10llllj) with ip^ (|21jl as input. The results are presented for two cases: when 
6 = 1 (the dominant term in asymptotics at |x| — > 00 is reproduced exactly, see ©) and when 6 
is taken as a variational parameter. In the latter the value of the second energy correction E2 is 
reduced in two times. 

the dominant and the first two subdominant terms in the asymptotic expansion © are reproduced 
by the function (|22H exactly. In this case the convergency condition (|15j) is satisfied, the data for 
Is/ilmaa: are i n Table II. These two cases we call the Case 1 (two parameters are fixed (see ([23))) to 
reproduce the growing terms in asymptotic expansion of y) and Case 2 (all five parameters in (|22jl 
are variational), respectively. Table II demonstrates that even the variational energy E^ 1 ' already 
provides extremely high accuracy: 10~ 8 - 10" 9 (Case 1) and 10~ 10 - lO" 11 (Case 2) for studied 
values of m 2 ,g. These results reproduce (or exceed) the best known numerical results in literature. 
The order of the third energy correction £^3 is already in the three-four orders of magnitude less 
than the second correction E<i- It indicates to the extremely fast convergence of the series in the 
parameter A (see ()11|0 , The energies E&)(= eW+E 2 + E 3 ) for given values of m 2 , g are the most 
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accurate among known in literature for the moment providing 16-17 significant digits (not all these 
digits are shown in Table II). 

In the Case 1 we are focused to reproduce the behavior of the wavefunction at large distances 
while the behavior at small distances is defined as a result of the variational procedure. We can 
calculate the coefficient in front of x in the logarithmic derivative of (|22|) which we denote E exp 
(see Table II). It should be the exact energy of the ground state in the case of the exact y (see 
0). It is worth emphasizing that the deviation of E exp from (or from the exact energy E) 
appears in the fifth significant digit (!). 

The important question concerns to the deviation of the function (1221) from the exact function. 
We consider the Case 1 when the growing asymptotics of y at large distances is reproduced exactly. 
Due to this fact the perturbation theory for any studied m 2 ,g is fast convergent and y\ is bounded 
for any real x, 

\Vl\max ~ 0.01 . 

Hence, the function (|22[) provides the uniform approximation of the exact y. In Fig.l for m 2 = 
— l,g = 2 the behavior of yo it is shown while Fig. 2 contains the behavior of the first correction y\. 
It is worth mentioning that the maximal value of y\\ \y-\\ m ax ~ 0.0029 appears at x ~ 3.9 where the 
value of ^0°'°^ * s already extremely small, ~ 10~ 9 and yo ~ 15.2. Then y\ oc 1/x 2 at i > 1. 
The ratio 2/1/2/0 ( see Fig. 3) is the extremely small function which tends to zero at |x| — > 00. A 
similar situation appears for other values of m 2 ,g. For the double-well potential while approaching 
the semiclassical limit: m 2 — > —00 and g is kept fixed, e.g. g = 1 - the accuracy provided by 
()22|) remains unchanged: E2 < —1. X 10~ 6 for m 2 ranging from -1 to -30. In the same range the 
parameters of (1221) are varied very little, being always of the same order of magnitude, although 
the value of the energy changes from ~ 1 to ~ —105. It was also calculated the parameters of the 
potential (1) which the ground state energy vanishes, 

E(m 2 = -2.2195970861,(? = 1) » 10~ 13 . 

It is worth mentioning that uniform approximation for eigenfunction of (2/c+p)th excited state, 

k = 0, 1, 2, . . . ,p = 0, 1 can be easily constructed as well. It has a form 

A + ax 2 /2 + bgx i /± \ 

{d 2 +gx 2 ) 1 / 2 J ' 1 ' 

where Pk is a polynomial of kth degree with positive roots. Parameters a, b are always chosen 
following (|23|) in order to reproduce the growing terms in the asymptotic expansion of y (Case 1 



x p P k (x z ) 

vT + c 2 x 2 



cxp 



b 

a = C j 6 + 771 


9 = 2 
m 2 = 1 


9 = 1 
m 2 = 


9 = 2 
m 2 = —1 


4/3 
1.79502293 
0.7796 


1.33813399 
- 
0. 


4/3 
0.37849775 
0.2535 


1.33361251 
- 
0. 


4/3 
1.10179879 
-1.8484 


1.33190626 

- 

-2.00747 


eW 


1.607541303542 


1.607541302751 


1.060362094762 


1.060362090514 


1.029560850845 


1.029560831475 


\yi j max 
Eexp Q 


0.0026 
1.607362918 




0.0029 
1.059963236 




0.0064 
1.030192509 




^3 
B (3) 


-0.107E-08 


-0.28E-09 


-0.428E-08 
1.060362090485 
-0.64E-12 


-0.30E-10 


-0.198E-07 
1.029560831059 
-0.47E-11 


-0.42E-09 


1.607541302469 


1.607541302469 


1.060362090484 


1.029560831054 


-0.12E-12 


-0.2E-13 


0.1E-14 


0.4E-13 


1.607541302469 


1.607541302469 


1.060362090484 


1.060362090484 


1.029560831054 


1.029560831054 














Ref.[4] 


1.607541303 


1.06036209 


1.029560848 



TABLE II: The Hamiltonian (1): Variational energy of the ground state and its corrections in perturbation theory (jlOHl 1|) with ip^ (j22j) as 

input (the values of the corrections £2,3 are rounded to very few first digits). For given m 2 ,g the first column corresponds to the parameters 

(2) 

in ^ which are chosen to reproduce exactly the first three terms in the asymptotic expansion ©, the second column all parameters are 
taken as variational. The energies with all accurate significant digits are framed. 
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FIG. 2: The first correction y\ as function of x for (|22|) 

for the ground state)). In order to find the other parameters the orthogonality constraint should be 
imposed: for fixed (k,p), (|24() should be orthogonal to previously constructed functions (m,p), at 
m = 0,1,2, ... ,k — 1. It fixes some parameters in (|24|) while the remaining parameters are used as 
variational. It is worth mentioning that very important physical characteristic of ^ at m 2 < is 
the energy gap between the energies of the ground state (0, and the first excited state (0, 1 in the 
semiclassical limit: m 2 is kept fixed and g — > oo. Although we do not provide concrete numerical 
results in this paper, this limit is reproduced with very high accuracy. 

As a conclusion I would like to recollect my first and, in fact, the only scientific encounter with 
F.A. Berezin. It was in mid-1970s just before his breakthrough results in supermathematics. I 
was a recent graduate of Moscow Institute for Physics and Technology just hired by Institute for 
Theoretical and Experimental Physics. Among theoretical physicists F.A. Berezin had a reputation 
of the unique mathematician who was able to understand what physicists are talking about and 
with whom one could discuss your mathematical difficulties. He had run a weekly seminar at 
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0.0002 




X 







-0.0002 







-0.0006 



-0.0010 



FIG. 3: The ratio y l /y for 



Mathematics Department of Moscow State University. By accident I went to one of his seminars 
but I could not say that I understood much. During the seminar break in the moment when I 
was ready to leave he approached to me, introduced himself and asked what I am working on. I 
answered that for already several years I had been trying desperately to find an analytic solution 
of the anharmonic oscillator x 4 (see Q). Or, at least, to understand why all my attempts had 
failed. He immediately said: "I know an anharmonic oscillator, which has an analytic solution! 
Take a function e~ x4 , differentiate it twice and divide the result by e~ xi . This is the anharmonic 
oscillator potential where you know an exact eigenfunction." Then we talked for several minutes. 
In the end, he said that together with M.A. Shubin (now a Distinguished Professor at Northeastern 
University in Boston) they would write a book on the Schroedinger equation and one day if I did 
not mind they could ask my opinion or advice. He introduced M.A. Shubin to me with whom I 
become a lifelong friend. After the seminar F.A. Berezin, M.A. Shubin and myself took a walk to 
the nearest metro station. I remember it was a very interesting conversation about mathematics 
but I do not remember what it was about - I was engrossed with thinking about the remarks of 
F.A. Berezin. 

This short meeting had struck me and influenced strongly my scientific life. From human point 
of view I was deeply impressed how serious and friendly he was towards me, almost still a student. 
After some time I understood that instead of solving the Schroedinger equation with an already 
given potential one can generate a zillion potentials for which a single eigenstate is known exactly. 
It led me to an idea of choosing one of such potentials as a zero approximation to construct a 
convergent iterative procedure for a given eigenstate, not for the whole spectra Another idea 
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was related to a natural question of how to construct a potential for which not one but two or 
more eigenstates can be found explicitly. It led to a discovery of quasi-exactly-solvable quantal 
problems J3|. A few years after our meeting F.A. Berezin tragically died and the chance to discuss 
with him disappeared, but all my life I keep a memory of him. Recently, I told this story to one 
of the greatest minds in today's mathematics who was among very close colleagues and friends of 
F.A. Berezin. He responded with sorrow that "... only after his death we realized how strong was 
his influence on all of us and, eventually, how great he was!" 
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